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Abstract 

In this paper, we present the next-to-leading order predictions for diphoton 
production in the ADD model, matched to the HERWIG parton shower using 
the MC@NLO formalism. A selection of the results is presented for d = 2 — 6 
extra dimensions, using generic cuts as well as analysis cuts mimicking the 
search strategies as pursued by the ATLAS and CMS experiments. 



1 Introduction 



Extra dimension models [TJ |2] are important new physics scenarios at the TeV scale 
that address the hierarchy problem and are being extensively studied at the LHC. 
Recently both ATLAS [3J and CMS [4J have looked for evidence of extra spatial 
dimension in the diphoton final state and put bounds on the fundamental Planck 
scale Ms in (4 + d)- dimensions, for the 7 TeV proton-proton collision. 

In the extra dimension models, Kaluza-Klein (KK) modes of a graviton (as a 
result of the graviton propagating in the extra dimensions) , could decay to a pair of 
photons. Interaction of the massive spin-2, KK modes hffi with the standard model 
(SM) particles localised on a 3-brane, is via the energy momentum tensor T^ v of the 
SM 
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where k = s/l&jr/Mp is the reduced Planck mass in 4-dimensions. In process involv- 
ing virtual exchange of KK modes between the SM particles, like in the diphoton 
production process, the sum of the KK mode propagator T>(s) is given by 

K *V{s) = « 2 V \— 
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where s is the partonic center of mass energy, d is the compactified extra spatial 
dimensions, A is the UV cutoff of the KK mode sum which is usually identified 
as Ms [6] and the integral I(A/Q) is given in [5]. Note that in Eq. El we have 
included the k 2 (suppression as a result of gravity coupling), which on summation 
over the high multiplicity of KK modes compensates the suppression in the ADD 
model. In this analysis we have followed the convention of [5]. 

Improved theoretical predictions to higher orders in QCD have now been per- 
formed for cross sections of pair production processes viz. di-lepton [TJ, di-gauge 
boson (77 [8], ZZ and W + W~ [9]), which in extra dimension models could result 
from the exchange of a virtual KK mode in addition to the usual SM contribution. 
The real emission of KK modes lead to large missing Et signals viz. mono-jet [TO] , 
mono-photon [UJ, mono-Z boson, and mono-iy ± boson [12]. The next to leading 
order (NLO) QCD corrections in some of the above processes are quite substantial 
and their inclusion in the computation also lead to a reduction of theoretical un- 
certainties, making it possible for the experiments to put more stringent bounds on 
the extra dimension model parameters. 

The diphoton final state is an important signal for extra dimension searches, as 
the branching ratio of a KK mode decay to diphoton is twice than that of a decay 
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to individual charged lepton pair. The quantitative impacts of the NLO QCD cor- 
rection to the diphoton final state for extra dimension searches have been studied 
in [8], where various infrared safe observable were studied using phase space slicing 
method. The factorisation scale dependence gets reduced when 0(a s ) corrections 
are included. Fixed order calculation truncated to NLO, at best yields results for 
sufficiently inclusive observable. Combining fixed order NLO and parton shower 
Monte Carlo (PS) [T3J Hlj . would extend the coverage of the kinematical region to 
consistently include resummation in the collinear limit and also make a more exclu- 
sive description of the final state and get as realistic as possible to the experimental 
situation. The flexibility to incorporate hadronisation models and capabilities to 
simulate realistic final state configurations that can undergo detector simulations 
are the main advantage for the experimental collaborations. 

ATLAS [3j and CMS [I] have analysed the diphoton invariant mass spectrum, 
using a constant K-factor for the full range of the invariant mass distribution to put 
lower bounds on extra dimension scale to NLO accuracy. However, this choice is not 
sensitive to possible distortions of distributions that can arise at NLO. Our present 
analysis will further help to put more stringent bounds on the model parameters. In 
this analysis we have considered various distributions for the ADD model parameters 
d — 2 to 6 with appropriate Ms value as bounded by the experiments Hj. In 
Table [U the Ms values for different extra dimensions d used in this analysis have 
been tabulated. For relevant observables we consider the fixed order results to NLO 
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Table 1: Ms values used for the various extra dimensions d in our analysis. 

accuracy and include PS. Factorisation, renormalisation scale uncertainties and PDF 
uncertainties are also estimated in an automated way [15]. For photon isolation, both 
smooth cone isolation and the experimental isolation criteria are considered. 

The rest of the paper is organised as follows. In section 2, we discuss the NLO 
results to the diphoton final state in the ADD model and the essential steps needed to 
implement the parton showering to NLO accuracy. In section 3, we present selected 
numerical results and estimate the various theoretical uncertainties. Finally, in 
section 4 our conclusions are presented. 

2 NLO + parton shower 

Since the KK modes couple universally to the SM particles through the energy 
momentum tensor, both the qq and gg channel would contribute to the diphoton 
final state at leading order (LO). In the SM, the gg channel starts only at next to next 
to leading order (NNLO) via the finite box contribution through quark loop and the 
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large gluon-gluon flux at the LHC makes this contribution potentially comparable to 
the LO results. In the invariant mass region of interest to extra dimension searches, 
the box diagram contribution is not significant enough [8]. 

All the partonic contributions to NLO in QCD have been calculated for the 
diphoton final state [8J, for both large (ADD) [1] and warped (RSI) [2] extra di- 
mension models. QCD radiative corrections through virtual one loop gluon and real 
emission of gluons to the gg->77 subprocess, would contribute to both SM and 
extra dimension models. The q(q) g — > q(q) 7 7 begins to contribute for both SM 
and extra dimension models at NLO. The LO g g — > 7 7 extra dimension process 
will also get one loop virtual gluon and real gluon emission radiative corrections. 
There will also be interference between the SM and extra dimension model to give 
contributions up to 0(a s ). In this analysis we have not included the 0(a s ) cor- 
rections as a result of the interference between the SM box diagram contribution 
and LO extra dimension contribution to the g g — > 7 7 subprocess, as it is quite 
suppressed in the region of interest to extra dimension models. 

The q(q) g — > q{q) 7 7 NLO contribution has an additional QED collinear sin- 
gularity when the photon gets collinear to the emitting quark and can be absorbed 
into the fragmentation function which gives the probability of a parton fragmenting 
into a photon. Parton fragmentation functions are additional non perturbative in- 
puts which are not very well known. At the LHC, secondary photons as a result of 
hadron decaying into collinear photons and jets faking as photon are taken care of 
by photon isolation criteria [3], H] which also substantially reduces the fragmentation 
contribution. Since the fragmentation is essentially a collinear effect, the fragmenta- 
tion function can be avoided by the smooth cone isolation proposed by Frixione |T2] , 
which ensures that in no region of the phase space the soft radiation is eliminated. 
The smooth cone isolation is able to eliminate the not so well known fragmentation 
contribution and at the same time, ensures infrared safe (IR) observable. Centered 
in the direction of the photon in the pseudo rapidity (rj) and azimuthal angle (0) 
plane, a cone of radius r = ^/Jrf— ?7 7 ) 2 + (0 — 7 ) 2 is defined. The hadronic activity 
H(r) is defined as the sum of hadronic transverse energy in a circle of radius r < r$ 
and Ej, is the transverse energy of the photon. For all cones with r <tq the isolation 
criterion H{r) < H(r) max has to be satisfied, where i?(r) max is defined as 



Efforts for the experimental implementation of the smooth cone isolation is on going. 

Automation is an essential ingredient of this work. We have chosen to work in 
the aMC@NLO framework [17], which automatises the MC@NLO formalism [13] to 
match NLO computations with parton showers. In this paper we present results 
matched to HERWIG [18]. For the NLO computation, isolation of IR poles and 
phase space integration are carried out by MadFKS [T9] . which automatises the 
FKS subtraction method [20] using the MadGraph [21] matrix-element generator, 
whereas for one-loop amplitudes the results of Ref. [8] are used. The automation 
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Figure 1: Transverse momentum distribution pip of the diphoton for the fixed order 
NLO and NLO+PS. The ADD model parameters used are d = 2 and Ms = 3.7 
TeV. The lower inset displays the scale and PDF, fractional uncertainties for the 
NLO+PS results. 



within the MadGraph framework requires a new HELAS [22J subroutine to calculate 
helicity amplitudes with massive spin-2 particles [231124] . In addition, for our present 
analysis we have implemented the sum over the KK modes to take care of the virtual 
KK mode sum (Eq. [2]) that contributes to process in the ADD model [21] . We use 
this framework to generate the events for 8 TeV run at the LHC For the invariant 
mass distributions we have reproduced the results of [8J using the fixed order results 
from MadFKS. Also numerical cancellation of the singularities from the real and 
virtual terms have been explicitly checked. 



3 Numerical results 

In this section, we present the results for various kinematic distributions of photon 
pair in SM and ADD model. We have included all the subprocess contributions to 
NLO. The following input parameters are used: = 132.507, Gp = 1.16639 x 10~ 5 
GeV" 2 , m z = 91.188 GeV and MSTW2008(n)lo68cl for (N)LO parton distribution 
functions (PDF) [25]. The MSTW PDF also sets the value of the strong coupling 
ots{ m z) at LO and NLO in QCD. The renormalisation and factorisation scales are 
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Figure 2: Invariant mass distribution M 77 for ATLAS (left panel) and CMS (right 
panel) for d = 2 extra dimensions and Ms = 3.7 TeV. The SM contribution to 
NLO+PS and ADD to LO+PS and NLO+PS have been plotted. For the NLO+PS 
(ADD) results, the lower inset displays the scale and PDF fractional uncertainties. 



chosen as \xp = lir = M 77 , the invariant mass of the photon pair. The events that 
have to be showered are generated using the following generation cuts: |^ 7l 2 | < 2.6, 
pip' 2 > 20 GeV, diphoton invariant mass 100 GeV < M 77 < Ms and the photon 
isolation is done using the Frixione isolation with r$ = 0.38, e 7 = 1 and n = 2 (see 
Eq. More specific analysis cuts can be applied subsequently to generate the 

events. 

The dependence of the prediction of an observable on the factorisation and renor- 
malisation scales, is a result of the uncalculated higher order contributions, which 
can be estimated by varying iip and /xr independently around the central value 
/J-f = /J-r = M 77 . The variation is done by the following assignment fxp = £p M 77 
and ix r = £r M 77 , where the values for (£f,£r) used are (1,1), (1/2,1/2), (1/2,1), 
(1,1/2), (1,2), (2,1), (2,2). The various ratios of fxp, lir and M 77 that appear as 
arguments of logarithms in the perturbative expansion to NLO are within the range 
[1/2,2]. The variation of both iip and llr are taken as the envelope of the above indi- 
vidual variations. Variation of only Lip would involve the choice £r = 1 and varying 
£f and vice-versa for variation of only lir. The PDF uncertainties are estimated 
using the prescription given by MSTW [25]. Fractional uncertainty defined as the 
ratio of the variation about the central value divided by the central value, is a good 
indicator of the scale and PDF uncertainties and is plotted in the lower insets to the 
various figures. As described in Ref. [15], the generation of these uncertainty bands 
can be done at virtually no extra CPU cost within the aMC@NLO framework. 

To begin with, we compare the fixed order NLO result with NLO+PS for the 
transverse momentum of the diphoton log 10 pip using generic cuts: M 77 > 140 GeV, 
|%| < 2.5, > 40 GeV, > 25 GeV and r = 0.4. In Fig.HJ log 10 pj 7 distribution 
is plotted for d = 2 with appropriate Ms value. It is clear that at low pip values, 
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Figure 3: Invariant mass distribution for d = 3 (left panel) and d = 4 (right panel) is 
plotted for ADD and SM contributions to NLO+PS accuracy. The lower insets gives 
the corresponding scale and PDF, fractional uncertainties for NLO+PS (ADD). 

NLO+PS correctly resums the Sudakov logarithms, leading to a suppression of the 
cross section, while the fixed order NLO results diverges for pjJ — > 0. At high p7p , 
the NLO fixed order and NLO+PS results are in agreement. In the lower inset of 
the Fig. [TJ we have the scale and PDF variation of the NLO+PS, which increase 
with as observed in [26] . 

We now present the results for the various kinematical distributions to NLO 
accuracy with PS (labelled as NLO+PS), for analysis specific cuts. Both the ex- 
periments ATLAS and CMS have looked for diphoton invariant mass in the region 
140 GeV < M 77 < M s . ATLAS cuts [3]: the rapidity of the individual photons 
are in the region |r/ 7 | < 2.37, with an exclusion region 1.37 < |r/ 7 | < 1.52, the 
transverse momentum of the individual photons p\ > 25 GeV and for photon iso- 
lation: sum of transverse energy of hadrons ^2-^t{H) < 5 GeV with Ar < 0.4. 
Ar = a/ A0 2 + At/ 2 is a cone in the azimuthal angle, rapidity plane. For CMS the 
corresponding cuts are jl]: |r/ 7 | < 1.44, p^ > 70 GeV, photon isolation: (i) sum of 
the energy of hadrons Yl-E(H) < 0.05-E 7 with Ar < 0.15, (ii) sum of transverse 
energy of hadrons Y. e t(H) < 2.2 GeV + 0.0025 with 0.15 < Ar < 0.4. In 
addition to the ATLAS and CMS photon isolation, if we also include the Frixione 
isolation criteria, there is no appreciable change in the results. 

In Fig. [21 we have plotted invariant mass distributions da/dM^ of photon pair 
in the SM as well as in the ADD model for ATLAS (left panel) and CMS (right 
panel). For ADD model we have obtained the distributions for Ms = 3.7 TeV and 
d = 2. The central value curves correspond to the choice Hf = Hr = M 77 , have been 
plotted for the ADD (NLO+PS) and purely SM (NLO+PS) contribution. The label 
ADD refers to the total contribution coming from SM, ADD and the interference 
between them. The corresponding ADD (LO+PS) contribution gives an indication 
of the quantitative impact of the NLO QCD correction. At larger invariant mass of 
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Figure 4: For the invariant mass distribution, with d = 2 and Ms = 3.7 TeV, the 
fractional uncertainties as a result of fip variation (left panel), [Ir variation (central 
panel) and //p, /ir variation (right panel). 




Figure 5: The rapidity (Y) distribution of the diphoton pair for d = 3 (left panel) 
and d = 4 (right panel) for SM (NLO+PS) and ADD (LO+PS and NLO+PS). The 
lower inset displays the corresponding fractional scale and PDF uncertainties of the 
NLO+PS (ADD) results. 

the photon pair the ADD effect is dominant. To demonstrate the sensitivity of our 
predictions to the choice of scale and PDF uncertainties, in the lower inset fractional 
uncertainty by varying (a) both fip and hr and (b) PDF error sets, are plotted. The 
difference in the distribution in Fig. |2] for ATLAS and CMS can be attributed to 
the very different cuts used for their analysis. In Fig. [3j the corresponding plots 
for d = 3, 4 are plotted for the CMS cuts. The choice of M s used for the plots 
corresponds to the lower bounds obtained by [3], H] using the diphoton process. 
By including higher order corrections, the scale dependence goes down from about 
25% at LO, to about 10% at NLO, as can be seen from the ratio plots. The PDF 
uncertainty does not change significantly and remains about 8%. 

We now consider the fractional scale uncertainties on the invariant mass dis- 
tribution as a result of the variation of the scales \ip and \ir (both independently 
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Figure 6: The transverse momentum distribution p^P of the diphoton for d = 3 (left 
panel) and d = 4 (right panel). 



and simultaneously) in going from LO+PS to NLO+PS. Note that the LO cross 
sections depend only on fxp through the PDF sets, but at NLO level the scale /Xr 
enters through g: s (/zr) and log(//p / '/xr) coming from the partonic cross sections af- 
ter mass factorisation. As expected the inclusion of NLO QCD correction reduces 
the factorisation scale dependence resulting from the LO observable which is clear 
from Fig. H] (left panel). In the high M 77 region, the uncertainty of about 25% at 
LO+PS gets reduced to 5% when NLO+PS corrections are included. On the other 
hand, the /xr dependence enters only at NLO level (see middle panel of Fig. Hj) 
which will get reduced only if NNLO corrections are included. Hence, we see our 
NLO corrections are sensitive to the choice of /i# but the variation is only 5% and 
is fairly constant for the range of invariant mass considered. If we vary both \ip and 
\Xr simultaneously as shown in Fig. H] (right panel), we find that the reduction in 
the ix p scale dependence at NLO level is mildly affected by the Hr variation in the 
large invariant mass region. In the small invariant mass region, the LO and NLO 
results exhibit smaller fxp dependence compared to the large invariant mass region. 
But lir dependence coming from the NLO results does not change much with the 
invariant mass M 77 . Hence variation due to lir at small M 77 is larger compared to 
that resulting from fxp- This explains the behavior at small invariant mass regions 
where the NLO+PS variation is in excess of the LO+PS (see right panel of Fig. H]). 

The rapidity distribution of the diphoton pair is plotted in Fig. [5] for d = 3 
(left panel) and d — 4 (right panel). For this analysis we have chosen M 77 > 600 
GeV, the region where the effects of ADD model begins to dominate over the SM 
diphoton signal at NLO (see Fig. [3]). The scale and PDF uncertainties to NLO 
are displayed as insets at the bottom of each figure. The scale uncertainties are 
usually larger than the PDF uncertainties in the rapidity distribution except for the 
central rapidity region where they are comparable. For d = 3 the scale uncertainties 
are about 20% around the central rapidity region, which come down to about 10% 
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when NLO+PS corrections are included. The PDF uncertainties for LO+PS and 
NLO+PS are comparable. 

Finally, we plot the transverse momentum distribution in Fig. O for d = 3 (left 
panel) and d = 4 (right panel), for the SM and ADD model to NLO+PS accuracy, 
with M 77 > 600 GeV. The ADD results are also plotted for LO+PS. The scale and 
PDF uncertainties are displayed as insets at the bottom of the plots for NLO+PS 
(ADD). 

4 Conclusion 

In this analysis, we have presented the diphoton final state in the large extra dimen- 
sion model to NLO in QCD and matching to parton shower is implemented using the 
aMC@NLO framework. All the subprocesses that contribute to the diphoton final 
state from both the SM and ADD model are considered to NLO in QCD. This is the 
first time MC@NLO formalism has been used for a processes in the ADD model and 
we hope it would significantly help extra dimension searches at the LHC to constrain 
the ADD model parameters. Using a set of generic cuts we first demonstrated the 
importance of NLO+PS over the fixed order NLO computations, by considering the 
pP distribution. We have presented our results for various observables viz., invari- 
ant mass, rapidity and transverse momentum of the diphoton, both for the ATLAS 
and CMS detector specific cuts to NLO+PS accuracy. It is important to note that 
there is substantial enhancement of the various distributions due to the inclusion 
of NLO corrections and both the theoretical and PDF uncertainties have been esti- 
mated. There is a significant decrease in theoretical uncertainties from over 20% at 
LO to about 10% when NLO corrections are included. The results are presented for 
different number of extra spatial dimensions d = 2 — 6 and the respective values of 
fundamental scale Ms that have been experimentally bounded. The event files for 
d = 2 — 6 are available on the website http:/ / amcatnlo.c ern.cn] and we are working 
on making the code that was used to generate these events publicly available. 
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